High resolution reflection shooting

ABSTRACT

An illustrative embodiment of the invention includes methods for obtaining seismic data from closely spaced subsurface reflection points using a minimal number of seismic energy generations. The invention includes methods of shot point and seismometer spacing which may be varied in a systematic manner to produce reflection from close, but assuredly not the same, subsurface bed interface areas.

Waited tates Patent Treyhig 1 1 Sept. 18, 1973 ENG RESOLUTION REFLECTION 3,286,228 11/1966 Anstey 340 155 CP SHOOTING [75] Inventor: James Alva Treybig, Houston, Tex. Primary EXami'IeF'BCn-iamiIP Borchelt Assistant Examinerl-l. A. Blrmrel [73] Assngnee: Texaco Inc., New York, N.Y. Atmmey ThomaS h l et [22] Filed: June 16, 1971 21 App]. No.: 153,554 [571 ABSTRAW An illustrative embodiment of the invention includes 52 U.S. c1 340/155 MC 340/155 CP meflwds 0bmining Seismic data fmm clsely SPacecl [51] int. Cl G0lv 1/20 Subsurface reflection Points using a mimmal number of 58 Field of Search 340/155 CP 15.5 MC Seismic energy genemms- The invention includes methods of shot point and seismometer spacing which [56] References Cited may be varied in a systematic manner to produce reflection from close, but assuredly not the same, subsur- UNITED STATES PATENTS face bed interface areas. 2,906,363 9/1959 Clay, Jr 340/155 CP 3,018,838 1/1962 Hammond 340/155 CP 13 Claims, 11 Drawing Figures SYMMETRY 3D 0 0 LINE D 1 *1 7 7" 1 F t RECORDER 5P4 (s/=1.) SP3 LISP?! SP2 5 s/=2) 5P1 F39 L i 37 38 3 4 46 33 32' \44 G3 G4 SHOT POINT GEOMETRY FOR HIGH RESOLUTION BED INTERFACE PATENTED SEP! 8 I975 SHEEI 2 F 3 FOLD 4FOLD CDP JAE

HwP R W EEK FIG. 4

HIGH RESOLUTION RECORD F/G. 5 D *1 i EARTH '3 1 SURFACE PATENTEU SEP1 8 I975 sum 3 or 3 FIG. 6B

m mm2omwmm TRACES TRA CES FIG.6E

TRA CES llillGll RESOLUTION REFLECTION SHOOTING This invention relates to seismic exploration methods and more particularly to methods of seismic exploration in which separate seismic energy reflection records are obtained from a plurality of closely spaced but not coincident reflection points on subsurface bed interfaces which data are then processed in a manner which amplifies desired primary reflections and which minimizes vibrations on the record data caused byother seismic phenomena.

Seismic exploration relates to methods of obtaining information regarding subterranean earth formations by transmitting mechanical vibratory or acoustic energy from a point at or near the surface of the earth downwardly into subsurface formations. Seismometers or geophones spaced apart on the surface of the earth generally along a line or other predetermined curve, sense reflected or refracted seismic energy and produce an output, generally in the form of electrical signals, which are recorded in analog or digital form, on magnetic tape for example, for further analysis.

A vibrating mass or an explosive charge located on the surface or in a shallow shot hole is typically used as a source of energy for seismic exploration. Vibratory or acoustic energy from this source travels outwardly and downwardly into the earth formations and is reflected from media interfaces or discontinuities beneath .the surface of the earth. The media interfaces or discontinuities causing energy reflection occur between media having acoustic impedance characteristics differing from each other. The impedance discontinuity causes a portion or component of the energy to be reflected toward the surface of the earth at each such interface. The reflected vibratory energies incident upon the geophones or seismometers at the surface are then converted into representative electrical signals which are amplified and recorded generally as a function of time and surface location of the seismometer. By timing the arrival of selected portions of the returning vibratory energy, valuable information can often be obtained regarding the depth, slope and other characteristics of subterranean earth formations.

Unfortunately, however, other vibrations or vibratory energy are normally present at the seismometer location which tend to obscure the recognition of desired reflected vibratory energy. For example, if a vibrating mass energy source is used at the surface, horizontally traveling acoustic surface waves in the upper crust layers of the earth are propagated. These surface Rayleigh or shear waves may have a component of their energy in the vertical direction. The seismometers, which are designed to discriminate against vibratory energy arriving from other than a vertical direction, are unable to distinguish the vertical component of the horizontally traveling shear waves from the vertical component of a compressional wave traveling generally in a vertical direction after being reflected from a subterranean earth formation. Moreover, these horizontally traveling waves may be reflected from surface media dicontinuities such as valleys or hills. Thus multiple reflections of surface waves resulting in standing waves or interference patterns can occur and be superimposed upon the record produced by the seismometer.

Similarly, portions of the compressional wave acoustic energy coupled from the energy source into earth formations near the surface may be refracted or bent at media interfaces in such a manner that they arrive at geophones spaced horizontally along the surface from the energy source at an angle such that they have a vertical component of energy which is sensed by the seismometer array. These arriving wave fronts, having been refracted, appear no differently to the seismometer from the desired reflected energy impinging vertically on the seismometer from deeper earth formations whose characteristics it is desired to record. Thus it is apparent that there are numerous sources of noise or undesired vibratory energy capable of introducing misinformation into the record produced by a given seismometer.

In the prior art, because of these and other sources of noise in seismic data records, it has become common I practice to utilize the common depth point or common reflection point method of obtaining seismic data as disclosed in US. Pat. No. 2,732,906. In this technique seismic energy sources are placed symmetrically about a center line between the energy source and the seismometer group in a manner such that, in the case of plane parallel subsurface interfaces, multiple reflections from substantially the same point or same small area on a media interface of subsurface earth formations are detected for each of a plurality of seismic disturbances created in performing the exploration. The multiple or redundant reflections from substantially the same reflection point impinging on the same geophone are then, after certain corrections, combined, as for example by summing, to produce a single seismic record from a given seismometer. In this record the seismic energy from a source has impinged upon substantially the same reflection point after traveling a plurality of different acoustic paths in the subterranean earth media. By averaging the data over these substantially different acoustical paths, but reflected from substantially the same point in the subterranean media, it is hoped to suppress random noise fluctuations while retaining true information indicative of the characteristics of media interfaces in the vicinity of the common subterranean reflection point.

In order to avoid the use of an undue number of seismic vibratory disturbances, multiple common depth point shots are generally made at substantially spaced intervals along the earth's surface. This introduces a shortcoming into the method, for example, in the case of discriminating against horizontally traveling vibratory energy having a vertical component such as the previously mentioned shear or Rayleigh waves. The common depth points occur at the geometric midpoint between the seismic shot point and the seismometer arrays (for the plane parallel layer case). If the common depth points are substantially spaced apart on the surface of the subterranean earth formation being investigated, it is effectively impossible to discriminate against or resolve surface disturbances having wave lengths of the order of the distance of the spacing between the common reflection points from each other on the subterranean formation. This results because of the phenomenon of aliasing or foldback.

The aliasing phenomenon arises from the fact that it is necessary to sample atraveling wave front at points occuring within its wave length in order to discern information concerning the shape of an individual wave cycle. Therefore, if surface waves travel along the upper layers of the earth between the seismometer groupings and have a wave length equal to or greater than the one-half the seismometer spacing, these waves may not be effectively cancelled by averaging the seismometer data in the manner commonly used in the common reflection point or common depth point technique due to the fact that an inconsistent sampling of the data is taken. Moreover, surface waves which have wave lengths comparable with the spacing of the seismometers may actually be reinforced or enhanced by the simple averaging technique used in the common depth point method because such waves can conceivably arrive in phase after each shot at a given seismometer location.

It is thus desirable to obtain information concerning subterranean earth formations in a manner which is not subject to the limitations of the common depth point method and which, moreover, may be used to obtain other and further information than is available by the use of the common depth point method. This may be accomplished using the principles of the present invention in which a plurality of longitudinally spaced geophone groups are used. According to the invention the shot point to geophone distance is varied unsymmetrically about the shot point/geophone center line during a series of systematically varied longitudinally spaced shots. From these shots (in the case of the parallel plane geometry) a series of close, uniformly spaced reflections from points along the acoustic interfaces of the subsurface formations are obtained. These closely spaced sample points on the subsurface interface are obtained in the present invention by spacing the shot points so that it is assured that a common depth point is not intercepted by the ray paths of theacoustic energy used to record the seismic data.

To obtain this effect in the present invention a systematic variation may be introduced in the shot point spacing. This shot point spacing variation insures reflection from points or small areas more or less uniformly spaced on the subsurface formation interface. In using roughly the same number of shot points per mile for generating acoustic energy the method of the present invention produces many times as many different sample points along the subsurface acoustic interfaces than are obtained with prior art techniques. This may be referred to as a high resolution data gathering technique.

In addition to obtaining more densely located data points of the subsurface earth formation reflection interfaces, the present invention includes improved techniques for enhancing the seismic data. These techniques include the application of spatial filters to seismic data obtained from the closely spaced subsurface reflection points. Spatial filtering techniques of varying degrees of complexity, according to the invention, are applied to seismic records or traces from adjacent densely spaced reflection points. This, in effect, averages in a weighted manner, data from a region in the vicinity of a given subsurface reflection point thus enhancing the validity of data at each such point.

The present invention, including other and further objects, features and advantages, is best understood by the following detailed description thereof when taken together with the accompanying drawings in which:

BRIEF DESCRIPTION OF THE DRAWINGS FIG. I is a schematic view showing the shot point geometry of the prior art common reflection point seismic method.

FIG. 2 is a schematic view showing one example of shot point geometry in the practice of the present invention.

FIG. 3 is a schematic view showing data enhancement techniques according to the prior art.

FIG. 4 is a schematic view illustrating data enchancement according to the present invention.

FIG. 5 is an illustrative view showing the effect of data gathering point geometry on signal aliasing.

FIGS. 6A, 6C, and 6E are graphical illustrations of spatial filter operators useable in practicing the present invention, and

FIGS. 68, 6D, and 6F are graphical illustrations of the effect on data of the operators of FIGS. 6A, 6C, and 6E.

DESCRIPTION OF THE PREFERRED EMBODIMENT Referring initially to FIG. 1, the present invention may be explained by reference to the shot point geometry of the prior art seismic exploration method known as common refection point or common depth point shooting. FIG. 1 illustrates this prior art shot point geometry. For the sake of simplicity it is assumed in the illustrations of FIG. 1, and as will be discussed subsequently FIG. 2, that only one subsurface earth formation interface exists and'that the surface of the earth is flat and parallel to the bed interface of the subsurface earth formation being investigated. In practice it will be realized, of course, that the subsurface geometry may be many times more complex than this.

In the prior art common reflection point method, a plurality of shot point locations, 11, 12, 13 and 14, for example, are located symmetrically about a line of symmetry 26. Symmetry line 26 falls between the shot points and the geophone or seismometer locations 15, 16, 17, and 18 for the case of flat, parallel plane surfaces as illustrated in the simplified model of FIG. 1. In practice it will be appreciated that while only one phone is shown at the various geophone locations, that a plurality of geophones can be used (in a geometric array if desired) to reduce the effects of surface noise. The principal aim of the common reflection point method is to reflect seismic energy from the different shot point locations off essentially the same points or same small areas (21, 22, 23, 24, 25 of FIG. l, for example) of the subsurface bed interface order investigation. The energy reflected from these points is then detected or sensed by the seismometers 6,, G G G etc. situated at the locations I5, 16, 17 and 18.

For ease of illustration, only four shot point and geophone locations are shown. It will be known to those skilled in the art that in practice 6, 12, or 24, or in fact any number of geophones and shot points may be used in common reflection point shooting, depending on the degree of accuracy or redundancy desired. If it is assumed that the velocity of propogation of acoustic energy in the media 20 is continuous over the spatial interval encompassing all shot points and geophone locations and if the geometry is as shown in FIG. I, then the symmetrical location of the shot points and geophones about the line of symmetry 26 results in principal paths (or so-called ray paths) of acoustic energy as illustrated. For example, if the seismic energy from shot point 11 traverses the ray path 11, 21, 15, it is sensed at the seismometer or detector 6,. Similarly, seismic energy from shot point 2, which is tired at a different time from the firing at shot point ll, traverses a path to geophone G (at location 16) over a ray path i2, 21, 116. This conclusion is reached as a result of applying the laws of reflection of acoustic energy at the interface between the subsurface beds. Moreover, energy from earlier or later shots located at shot points 3 and 4t (following the paths 13, 2T, 17 for the case of shot point 3 and 1d, 21, lid for the case of shot point 1) are all reflected from point 2f at the bed interface of the subsurface formation. Thus, four shots at shot points ll, 2, 3, and 4 of FlG. il produce four reflections which are received at four geophones G G G and G all of which had been reflected from a common depth point or the same small area of the bed interface 21. Recorder 113 records signals from the geophones as a function of time. lily repeating shots at other analogously placed shot points and receiving reflections at analogously placed geophone locations, fourfold redundant data may be gathered from each reflection point.

The common reflection point shooting technique then assumes that the data reflected from the same or common depth point, as for example, point Zll, is in some sense related and thus may be algebraically combined, as for example by addition, to yield an enhanced signal to noise ratio for the true seismic energy reflected from the subsurface bed interface being examined. it is further assumed that noise, being of a random (or at least less consistent) nature, will be disenhanced or dispersed by the averaging or addition technique.

It will be appreciated by those skilled in the art that before the averaging technique is applied that the seismic data must be corrected for moveout. Moveout is the apparent time shift of seismic events at horizontally spaced geophone locations caused by the horizontal spacing of the geophones and the finite velocity of the acoustic energy. There are many schemes known in the art for doing this. Moveout methods will not be discussed herein. It will suffice to say that whenever the averaging or spatial filtering technique is discussed it is assumed that the data has been corrected for moveout.

As previously discussed, ray paths for seismic energy detected at the seismometers G G G G may be traced for various combinations of shot points and geophone or seismometer locations, as for example the other paths shown in FIG. ll. The resultant common depth points or common reflection points on the subsurface bed interface being investigated (if the geophones-G G G and G, are located respectively a distance D apart) are a series of subsurface reflection points 211, 22, 23, 24, 25 etc. which are located a distance D/2 apart.

it will be noted that using-the common depth point or common reflection point scheme of gathering seismic data that the reflection points on the bed interface in question are generally a relatively long distance apart. This results from the practical consideration that in order to cover ground the shot points and geophones are located a relatively long distance D apart, typically, 200 to 300 feet. In cases of practical interest this can result in the undesired gathering of noise data due to the phenomenon of aliasing because of the relatively sparse locations of the data points on the subsurface bed interface being investigated.

Referring now to FIG. 5, the phenomenon of spatial aliasing is illustrated schematically. it is well known that a surface wave will be generated in the layer of earths media nearest the surface from the release of the energy for producing the seismic disturbance. This surface wave will have components of energy in the vertical direction which may be sensed, for example, by the seismometers or geophones located at locations 115, lb, f7 and T3 of FIG. ll. Points 15, lib, l7, and T8 in H6. 5 correspond to the same numbered points in FIG. ll. if the frequency content of the surface wave is sufficiently high, it may have an overall wave shape such as that shown by the solid line wave in FIG. 5 of wave length it. The geophone locations, on the other hand, are in this case located more than one wave length A of the surface wave apart, i.e., here D it Therefore, for any given instant of time the geophones located at 15, llti, 117 and lid may sample amplitude characteristics of the surface wave in the manner shown in FIG. 5 by the circled points. In this example, the seismometer at location 15 may sample a slightly negative amplitude characteristic value. The seismometer at location 116 may sample a slightly more negative value than that at 115. That at location i? may sample nearly the same negative value as that at location 115. Finally, the seismometer at location 113 may sample a slightly positive value. This data sampling would seemingly imply that a surface wave having a wave length similar to that shown by the dotted line curve in FIG. 5 could be present. It is thus apparent that because of the spacing distance D of the geophones being greater than the wave length It of the traveing surface wave that a false impression has been created. This effect is known as aliasing. Aliasing can result in apparent wave length stretching, phase shift, wave fold back or other apparent distortion of the data.

It is apparent from the example of FIG. 5 that whenever the sampling distance D is greater than one-half the wave length A of the traveling surface disturbance, for example, that the possibility exists that the surface disturbance will be incorrectly interpreted due to the aliasing phenomena. With the use of the present invention it is possible to achieve, using the same number of seismic shots or energy generations, a much denser or closely spaced field of investigation on the subsurface bed being investigated than is possibleusing the common reflection point or common depth point technique. This can be used to advantage in overcoming aliasing as will be described.

Referring now to FlG. 2, a simplified sketch of the shot point geometry in accordance with the present invention is shown. In the example of FIG. 2, again it is assumed for simplicity that the surface of the earth is flat and that the bed interface being investigated is also flat and parallel to the earths surface. Again geophones or seismometers G G G G, are located at locations 35, 36, 3'7 and 3b., spaced equally a distance D apart. The seismometer G is spaced a distance D from the line of symmetry 39. The initial shot point in a sequence of shots is at location 3ll. This is located a distance D to the left of the line of symmetry 39. However, in the present invention, shot point 2, instead of being located a distance D to the left of shot point ll (i.e., at location Ml of FIG. 2 which would be its common depth point location) is instead located at point 311 a distance D D/4 to the left of shot point ll. Shot point 3 is located a distance D D/d to the left of the location of shot point 2 (32 of HG. 2) and thus falls at location 33 of FIG. 2 rather than at location 45 (which would be the location of shot point 3 in the common reflection point technique.) Finally, shot point 4, which completes a shot sequence analogous to the fourfold redundant common reflection point technique of FIG. I, is located a distance D 3D/4 to the left of the psuedo (common depth point) location point 45 of shot point 3. Thus in the technique of the present invention (in the case of fourfold redundancy shooting corresponding to the fourfold redundance common reflection point shooting illustrated in FIG. ll) the sequence for shot point locations is D/4, 2D/4, 3D/4, to the left of each analogous common depth point shot location. For n-fold redundant shooting the locations of shot point would be D/n, 2D/n, (n l) D/n from the corresponding analogous n-fold CDP shot location.

The result of locating the shot points in this manner is to produce a relatively closely spaced series of reflection points on the subsurface bed interface being investigated. For example, the seismic enegy generated by the first shot at shot point l and detected by geophone G at location 35 follows the ray path 31, 451, 35. A reflection point is established on the subsurface bed interface at 41. Energy from the same shot, reaching geophone G at location 37 traverses a ray path 3K, 43, 37. Energy from the second shot (located at 32), rather than traversing the common depth point all in reaching geophone G at location 36, traverses the path 32, d2, 36. Here point 42 is located a distance D/8 to the left of point (it. This is due to the displacement of shot point 2 (located at 32) from its pseudo (or corresponding CDP location, 44?). Using this spacing technique it may be shown that a plurality of reflection points each spaced a uniform distance D/8, for example, on the subsurface interface may be produced. It will be understood that this technique may be used to obtain even more densely spaced reflection points. For example, if it is desired to use spacing analogous to sixfold CDP shooting, then using the method of the present invention spacings would be D D/6 apart (for those after the first shot) and would yield reflection points D/l2 apart on the subsurface interface being investigated.

According to the principles of the present invention, the data gathered from these finitely spaced reflection points on the subsurface bed interface are then combined by improved techniques to produce an enhanced data record. It is apparent from the geometry of FIG. 2 that if all of the ray paths are traced from the various shot points to each of the geophones in FIG. 2 a succession of relfection points spaced a distance D/S apart on the subsurface bed interface are produced. The energy from each of the shots has traversed a path through each of these points on its way to the geophones and no energy has passed through the same points or substantially the same points as in the case of the common reflection point technique. These systematically spaced reflection points are, however, spatially related to each other on the subsurface bed interface so as to be combinable according to techniques to be subsequently described. Such combination or processing can enhance the signal to noise ratio of the data in a manner analogous to that of the common reflection point technique while at the same time providing a high resolution, closely spaced image of the subsurface bed interface being investigated.

Referring now to FIGS. 3 and 4, the combination of the data to enhance its signal to noise ratio as done in the prior art common depth point technique and with the method of the present invention are contrasted. In FIG. 3 which could correspond to the example shown in FIG. ll, four separate seismic traces resulting from energy passing through each common reflection point on the subsurface formation are available. These are illustrated schematically by the four curves 51, 52, 53, and 54, for example, on FIG. 3. These data curves which have resulted from energy passing through the same point are then combined by algebraically adding the four curves together. It is assumed that noise will appear as more or less random variations on each of the four traces and will then cancel while the true signal resulting from the reflection of seismic energy at the common depth point will be systematic in nature and will be enhanced when the four traces are added together. Thus a resulting trace 55 is obtained in which the reflection from the subsurface bed indicated by the peak 56 on the curve is much enhanced and stands out over the surrounding noise on the trace.

Now referring to FIG. 4, a similar and yet much more powerful function may be performed to enhance the signal to noise ratio of the data gathered from the shot point geophone location procedure of the present invention. In this case the record of FIG. 4 may be thought of as a multiplicity of traces 61, 62, 63, 64, 65, 66, 67 etc. each generated by a separate but systematically spaced reflection point or data point located on the subsurface bed interface in question. Assuming that the moveout correction has been made, then according to the present invention the resultant signals 61, 62, etc. may be combined by means of a spatial filter in order to enhance the signal to noise ratio.

It will be appreciated by those skilled in the art that while FIG. 2 has illustrated a method for achieving a systematic variation in the location of subsurface reflection points or data points by monotonically increasing the shot point spacing, that using the data enhancement techniques to be described, this systematic variation of shot points is not required. It may, for example because of terrain, be impossible or difficult to place shot points in this manner. In this case the data enhancement techniques to be described would be equally applicable to non-uniform or nonsystematically spaced relfection point data provided the relative shot point/geophone location is such as to provide reasonable proximity of the reflection points or data points to each other on the subsurface interface being investigated.

The spatial filtering technique of the present invention is a means of averaging over the spatial dimension data resulting from reflections of acoustic energy from associated points on the subsurface bed interface being investigated. If, as is assumed in the common reflection point technique, the data resulting from reflected energy arriving from adjacent longitudinally spaced subsurface points to geophones (after being corrected for moveout) is in some sense associated, and if the noise source is more random in nature than the data source, then the signal to noise ratio of such data may be enhanced by an averaging technique. In the case of the common reflection point data, as previously discussed, the averaging technique is simply addition of the signals.

Referring now to FIG. 6A, a spatial filter operator is shown which corresponds to the averaging operator of the common depth point or common reflection point seismic method. The vertical scale of FIG. 6A, is taken as being indicative of the weight for effect of the data applied to the average by each of the data points. The

horizontal scale represents the location of the reflection point or data point in the spatial dimension. The filter operator expressed in the space domain is applied to the data resulting from energy reflections arriving from the spaced points labeled 2, -l, 0, +1, and +2 in FIG. 6A by taking the value of the data curve at the time and space point in question, multiplying it by the weight function bed interface being investigated are in general located much closer together than points obtained by the common depth point technique. This permits the removal of noise in a more effective manner than previously possible because the effect of spatial aliasing is much diminished. Moreover, more elaborate spatial filtering techniques may be applied to the data if desired.

The type of space domain operator illustrated in FIG. 6A is a linear operator in the sense that contributions from an arbitrarily wide region (determined byL, the width of the filter operator) are assumed to equally effect the filtered or resultant trace which occurs at the center of the operator and to contribute equal amounts to the resultant trace produced by the operator. If the effect of this type of linear averaging on data of a sinusoidal nature is examined mathematically, 'a resultant curve such as that shown in FIG. 6B is obtained. The curve of FIG. 6B is obtained by taking the Fourier transform of the space domain operator of FIG. 6A. That is, if the space domain operator is represented by f(x), where x is the spatial dimension along the line of the geophone and shot point locations, then the Fourier transform F(w) of this space domain function f(x) is given by:

F(W) a dx=fff ma dw (1) where W (Zn/A) and A is the wave length of the sinusoidal wave.

By taking the Fourier transform of the space domain operator and observing it in the wave length domain, the effect of the space domain operator on sinusoidal type functions of various wave lengths may be easily observed. (It may be assumed that most seismic data will fall into this class of functions.) The percentage of contribution of a given point to the resultant of the space domain operator is plotted as a function of the dimensionless quantity L/lt as in FIG. 63. It is seen that (for a fixed width L of operator) that as long as the wave length of the sinusoid being averaged is long (i.e., Ll) small) then the averaging technique is successful. However, as the ratio L/lt increases (i.e., as the-wave length of the data decreases), for example, it is seen that considerable distortion of data can occur when the linear average operator is used. Going further to the right on the curve of FIG. 68, it will be observed that for L/). ratios between It and 3/2, a 180 phase shift actually occurs by the application of the space domain filter of FIG. 6A. This is the foldback or aliasing effect previously described. If L is made larger (to better enhance the data by taking more points into account) then aliasing can distort the data even more with the linear average type of spatial filter operator. From this discussion it is apparent that some better form of spatial filter operator than the linear average operator of FIG. 6A, and the prior art CDP method, would be desirable.

. 1. That at 3 traces to the left of the 0 trace receives a weight of 2, etc. The central trace itself receives a weight of 5. The operator of FIG. 6C could operate on the high resolution record data of FIG. 4, for example, simultaneously on 9 traces corresponding to 9 equidis tant spaced points on the subsurface bed interface being investigated. In FIG. 4, for example, these could be traces 61 through 69. Again the resultant of the operator is a trace which may be ascribed to the central location (or 0 point) of the operator. The operator 0 point location is moved and the entire operation repeated to obtain continuously enhanced data.

Examining the response of the space domain operator of FIG. 6C (when operating on sinusoidal type functions), the Fourier transform resultant is plotted in the wave length domain in FIG. 6]). It is seen that the operator of FIG. 6C when applied to the data has the effect of substantially smoothing the response as compared to that of the operator of FIG. 6A. In this case it is seen, for example, that for a given L, a shorter wave lengths may be used in the average before distortion occurs. Moreover, in the caseof the operator of FIG. 6C, it is seen that no phase reversal occurs at any wave length as is possible in the case of thelinear operator of FIG. 6A.

An even more effective space domain operator is illustrated inFIG. 6E. In this case the operator length 'L has a width of 17 traces as illustrated by the trace numbers 8 to +8 on the horizontal axis of FIG. 6E. The weighting function applied to each of the traces in question is illustrated by the vertical scale in terms of a percentage. Here the (I point or central point of the operator corresponds to a weight of while traces removed as far as 8 space point locations from the center on either side of the center trace are attributed a weight of only I in summing the traces to produce the resultant. The response of this operator in the wave length domain is seen in FIG. 6F. Virtually no distortion is introduced by the short wave lengths (Ll). I) in the use of this operator.

From the discussion of the space domain filters illustrated here, it will be apparent to those skilled in the art that other more complicated space domain filter operators may be used. For example it may be possible to design such operators to perform specific functions if desired. Analog computer circuitry may be used to perform the space domain filter operator functions discussed above. Techniques for implementing such circuitry are known in the seismic art. For example, once the space domain operator is chosen, the spatial filter functions are represented by addition and multiplication. Thus amplifiers may be used which have gains proportional to the weighting function components and analog summing circuits known in the art used to sum the output of a plurality of such amplifiers for performing the spatial filter function on analog field data. Alternatively if the 'data is digitally recorded in the field, a general purpose digital computer may be proii i grammed to operate on the data at a later time in order to perform the space domain filering. The general operative scheme of such a program is apparent from the FIGS. 6A, 6C, and 6E. It is important to recall that the space domain filtering or averaging techniques are applicable to the data gathered by the high resolution method of the present invention in a much more useful manner than that applied to the prior art common reflection point shooting data due to the fact that reflection points or data points occur at much closer intervals on the subsurface bed interface being investigated than do the common reflection points. This enables noise of much shorter wave lengths to be filtered from the data without undue interference from the aliasing effect previously discussed.

The examples shown in the drawings are submitted as being illustrative only of the concepts of the present invention. While only four spaced shot points have been illustrated for simplicity in FIGS. 1 and 2, in discussing the geometry of the common reflection point technique and the high resolution data gathering technique of the present invention, it will be appreciated by those skilled in the art that as many shot points and geophone locations as desired may be used. For example, it is standard practice in the common reflection point technique to use 6, 12 or even 24 fold redundancy in gathering data points from the common reflection points on the subsurface bed interface being investigated. Similarly, in. the present invention, rather than using fourfold separation of the shot points from their pseudo location based on the common reflection point technique, it is impossible to use n-fold separation to achieve evenmore closely spaced data reference points on the subsurface bed interface being investigated. Moreover, while shot points have been shown only to the left of geophone locations, it is well known in the art to space geophones on either side of a shot point, and, as the shot points progress to the left, for example, to move portions of the geophone cable from the right end to the left end to perform a continuous survey. Similarly, rather than moving shot points to attain the geometry relationships desired, the geophones may be moved while the shot point is fixed. These same techniques of continuous survey are, of course, applicable to the method of the present invention as long as the basic geometrical relationships for systematically varying the shot point location as illustrated in FlG. 2 are employed in gathering the data. Finally, the data enhancement techniques, particularly with weighted filter operator, may be applicable to data gathered using nonsystematically spaced shot points, if desired.

While a particular embodiment of the present invention has been shown and described, it will be appreciated by those skilled in the art that other equivalent techniques may be suggested. it is the aim in the appended claims to cover all such changes and modifications as fall within the true spirit and scope of the present invention.

What is claimed is:

l. A method for obtaining seismic data records in seismic surveying comprising the steps of:

locating a seismometer and a first source of seismic energy in longitudinally spaced relation to each other on or near the earth's surface;

generating seismic energy from said first source and receiving at the seismometer reflected seismic energy components resulting from said energy generation and generating first electrical signals representative thereof;

relocating the seismometer and a second source of seismic energy in a different longitudinally spaced relation relative to each other at points so located that seismic energy propagated by said second source will travel a path of different length but will be reflected from an intermediate small area on a sub-surface bed interface which is a distinctly separated distance from the corresponding intermediate small area on the same subsurface bed interface from which reflected seismic energy components from said first source of seismic energy were reflected; and

generating seismic energy from said second source and receiving at the seismometer reflected seismic energy components resulting from said energy generation, generating second electrical signals representative thereof and recording said first and second representative signals.

2. The method of claim 1 wherein the steps of relocating the seismometer and a second source of seismic energy in a different longituidnally spaced relation to each other are repeated with each relocation of the seismometer and subsequent source of seismic energy being located in a systematically increasing longitudinally spaced relation.

3. The method of claim 2 wherein the subsequent sources of seismic energy are located in monotonically increasing longitudinally spaced relation to each other thereby producing monotonically increasingly distant small intermediate areas for reflection of seismic energy components on subsurface bed interfaces.

4. A method for obtaining seismic data records in seismic surveying comprising the steps of:

locating a line of seismometers in longitudinally spaced relation to each other on or near the surface of the earth;

locating a plurality of shot points for generating seismic energy on or near the surface of the earth in a line and in systematically increasing longitudinal spaced relation to each other and to a given seismometer;

generating seismic energy at each of the shot point locations;

receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and

recording said representative signals.

5. The method of claim 4 wherein said shot points are located in monotonically increasing longitudinally spaced relation to each other.

6. The method of claim 4 wherein said seismometers are located a distance D apart and if n shot points are located in longitudinally spaced relation, each shot point spaced subsequent to a first shot point in the line is located a distance D D/n from the precedingly spaced shot point.

7. A method for obtaining seismic data in groups of n fold associated seismic data records comprising the steps of:

locating a plurality of seismometer locations in a line and in equidistant horizontally spaced relation D to each other on or near the surface of the earth; locating n shot point locations in a line and in systematically increasingly distant longitudinally spaced relationship to each other and to a given seismometer on or near the surface of the earth;

generating seismic energy at each of the shot point locations;

receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and

recording said representative signals as a function of time after the occurrence of each energy generation.

8. The method of claim 7 and further including the step of relocating said plurality of seismometers after n such energy generations and repeating the succeeding steps for each such relocation to obtain a plurality of such n fold associated seismic data records.

9. The method of claim 7 wherein each of the spaced shot point locations is located a monotonically increasing longitudinally spaced distance, D D/r z, from its precedingly spaced shot point in the group of n shot point locations. 7

10. A method of obtaining seismic data'records in seismic surveying comprising the steps of:

locating a line of seismometers in longitudinally spaced relation to each other on or near the surface of the earth;

locating a plurality of shot points for generating seismic energy on or near the surface of the earth in a line and relative to said seismometers such that the spacing distance between each shot point and the mid-point equidistant between the shot points and the seismometers or the distance between each seismometer and said mid-point increases in a predetermined systematic manner; generating seismic energy at each shot point location; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals. ill. The method of claim 10 wherein said spacing distance increases in a monotonically increasing manner. 12. A method for obtaining seismic data records in seismic surveying comprising the steps of:

locating a line of shot points for generating seismic energy in longitudinally spaced relation to each other on or near the surface of the earth; locating a line of seismometers in systematically increasing longitudinal spaced relation to each other and to a given shot point; generating seismic energy at each of the shot point locations; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals. 13. The method of claim 12 wherein said systematically increasing longitudinal spaced relation comprises a monotonically increasing spaced relation.

(5w UMTED STATES PATENT OFFICE CERTIFICATE OF CORRECTFON Patent No. 3,760,347 Dated September 18, 1973 Invented s) JAMES A. TREYBIG I It: is certifie'd'that error appears in theabove-identi-fied patent and that said Letters Patent are hereby corrected as shown below:

' Column 6-, line 28, the word "traveing" should be --traveling--.

Column 6, line 61, the number "31" should be --32-.

Column 7, line W, "relfection" should be -reflection--.

Column 9, line 8, after the word "function" insert the following paragraphs which were omitted: w

of the filter operator at that space point and summing this with the results of the same operation performed on the data from each of the other space points of the space domain operator. The resultant is then associated with the space location falling at location 0 of the operator of Figure 6A.

The operator is then shifted to the next space location to be operated on (i.e., this point is placed at location 0' of Figure 6A) and the entl re process is repeated. In applying the operator of Figure 6A, for example, to the traces of Figure 3 for GDP gathered data, (assuming that there are fivefold CDP data present on the record for convenience) the curves 51, 52, 53, 5 and 55 would correspond to the data taken at the trace operator's data points'labeled -2, -l, O, l, and 2. Each data value on a curve would be multiplied by the weighting factor (1,

in this case) and then summed together The resulting 53 3 UNJETED STATES PATENT OFFICE Page 2 CERTIIFICATE OF CORRECTFON atent N 3,7 'J LT Dated September 18, 197% Inventor(s) JAMES A. TREYBIG It is certified that error appears in the abov e-identified atent and that said Letters Patent are hereby corrected as shown below trace is then associated with the point located in the center of the space domain operator. It is evident that this process is equivalent to simply adding all the curves together as previously described in processing GDP data.

When applied to seismic data obtained by the method of the present invention and represented, for example,

by the traces of Figure A, the Figure 6A spatial filter operator would, for example, take data at corresponding 7 times from curves 61, 62, 63, 6 65 and sum these data together for all corresponding times. The resultant curve would thus represent spatially filtered data occurring at the space dimension location of the curve 63 of the data records It will be recalled that in the method of the present invention, the reflection points from the subsurface-- Column 11, line 2 "filering" should be --filtering--. Column 12, line 23, 'longituidnally" should be longitudinally--. I

Signed and sealed this 9th day of April 197L (SEAL) Attest:

EDWARD LLFLETGHERJR. G MARSHALL DANE Attesting Officer Commissioner of Patents 

1. A method for obtaining seismic data records in seismic surveying comprising the steps of: locating a seismometer and a first source of seismic energy in longitudinally spaced relation to each other on or near the earth''s surface; generating seismic energy from said first source and receiving at the seismometer reflected seismic energy components resulting from said energy generation and generating first electrical signals representative thereof; relocating the seismometer and a second source of seismic energy in a different longitudinally spaced relation relative to each other at points so located that seismic energy propagated by said second source will travel a path of different length but will be reflected from an intermediate small area on a subsurface bed interface which is a distinctly separated distance from the corresponding intermediate small area on the same subsurface bed interface from which reflected seismic energy components from said first source of seismic energy were reflected; and generating seismic energy from said second source and receiving at the seismometer reflected seismic energy components resulting from said energy generation, generating second electrical signals representative thereof and recording said first and second representative signals.
 2. The method of claim 1 wherein the steps of relocating the seismometer and a second source of seismic energy in a different longituidnally spaced relation to each other are repeated with each relocation of the seismometer and subsequent source of seismic energy being located in a systematically increasing longitudinally spaced relation.
 3. The method of claim 2 wherein the subsequent sources of seismic energy are located in monotonically increasing longitudinally spaced relation to each other thereby producing monotonically increasingly distant small intermediate areas for reflection of seismic energy components on subsurface bed interfaces.
 4. A method for obtaining seismic data records in seismic surveying comprising the steps of: locating a line of seismometers in longitudinally spaced relation to each other on or near the surface of the earth; locating a plurality of shot points for generating seismic energy on or near the surface of the earth in a line and in systematically increasing longitudinal spaced relation to each other and to a given seismometer; generating seismic energy at each of the shot point locations; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals.
 5. The method of claim 4 wherein said shot points are located in monotonically increasing longitudinally spaced relation to eaCh other.
 6. The method of claim 4 wherein said seismometers are located a distance D apart and if n shot points are located in longitudinally spaced relation, each shot point spaced subsequent to a first shot point in the line is located a distance D + D/n from the precedingly spaced shot point.
 7. A method for obtaining seismic data in groups of n fold associated seismic data records comprising the steps of: locating a plurality of seismometer locations in a line and in equidistant horizontally spaced relation D to each other on or near the surface of the earth; locating n shot point locations in a line and in systematically increasingly distant longitudinally spaced relationship to each other and to a given seismometer on or near the surface of the earth; generating seismic energy at each of the shot point locations; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals as a function of time after the occurrence of each energy generation.
 8. The method of claim 7 and further including the step of relocating said plurality of seismometers after n such energy generations and repeating the succeeding steps for each such relocation to obtain a plurality of such n fold associated seismic data records.
 9. The method of claim 7 wherein each of the spaced shot point locations is located a monotonically increasing longitudinally spaced distance, D + D/n, from its precedingly spaced shot point in the group of n shot point locations.
 10. A method of obtaining seismic data records in seismic surveying comprising the steps of: locating a line of seismometers in longitudinally spaced relation to each other on or near the surface of the earth; locating a plurality of shot points for generating seismic energy on or near the surface of the earth in a line and relative to said seismometers such that the spacing distance between each shot point and the mid-point equidistant between the shot points and the seismometers or the distance between each seismometer and said mid-point increases in a predetermined systematic manner; generating seismic energy at each shot point location; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals.
 11. The method of claim 10 wherein said spacing distance increases in a monotonically increasing manner.
 12. A method for obtaining seismic data records in seismic surveying comprising the steps of: locating a line of shot points for generating seismic energy in longitudinally spaced relation to each other on or near the surface of the earth; locating a line of seismometers in systematically increasing longitudinal spaced relation to each other and to a given shot point; generating seismic energy at each of the shot point locations; receiving reflected components of seismic energy resulting from each such energy generation at each seismometer location and generating electrical signals representative thereof; and recording said representative signals.
 13. The method of claim 12 wherein said systematically increasing longitudinal spaced relation comprises a monotonically increasing spaced relation. 